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Abstract. The focus of this study is on the treatment 
of those stellar regions immediately adjacent to convec¬ 
tive zones. The results of hydrodynamical simulations by 
Freytag et al. (1996) show that the motion of convective 
elements extends well beyond the boundary of the convec- 
tively unstable region. We have applied their parametrized 
description of the corresponding velocities to the treat¬ 
ment of overshoot in stellar evolution calculations up to 
the AGB (Pop.I, Mzams = 3Mq). 

Our calculations show the 3'''^ dredge-up already at 
the 7*** thermal pulse (TP), and the dredge-up parameter 
reaches A = 0.6 during the next five pulses. Accordingly, 
the amount of dredged up is up to 10“^Mq. Our 
models develope a small so-called ^^C-pocket consisting 
of a few Finally, this treatment of boundaries of 

convective regions leads to intershell abundances of typi¬ 
cally ('^He/^^C/^®0)=(23/50/25) (compared to (70/26/1) 
in the standard treatment). 
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1. Introduction 

The question of burning and mixing along the AGB is in¬ 
timately linked to the so-called carbon star mystery (Iben 
1981): Observations of the LMC and SMG prove that most 
of the carbon (and s-process element) enriched AGB stars 
are found at rather low luminosities indicating that they 
also have rather Zow masses (1...3Mq, cf. Smith et al. 1987, 
Frogel et al. 1990). On the other hand, AGB models pre¬ 
dict the 3“^'^ dredge-up, i.e. the mixing of interior carbon 
to the surface during the recurrent He shell flashes on 
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the AGB, mostly for much higher luminosities, i.e. larger 
(core) masses. Thus, we meet two problems: (i) Why do 
we not observe luminous carbon stars?; and (ii) Why do 
we not find dredge up in low mass AGB models? 

The first one belongs most likely to the occurrence of 
hot bottom burning and mass loss (Iben 1975, Blocker 
& Schonberner 1991, Boothroyd et al. 1993), see also 
D’Antona & Mazzitelli (1996) for a recent discussion. 

The second problem will be subject of this Letter. 
Synthetic AGB calculations utilize the stellar parame¬ 
ters known from evolutionary AGB model sequences and 
take the minimum core mass for dredge up, and 

the dredge-up parameter A (ratio of dredged-up mass 
to growth of the hydrogen exhausted core per flash cy¬ 
cle) as adjustable parameters (among others) in order to 
match the observations, i.e. the luminosity function of 
LMG carbon stars. For example, van den Hoek & Groe- 
newegen (1997) find Mg"" = O.58M0 and A = 0.75 as 
best fit in contrast to the results of evolutionary calcula¬ 
tions (Mf'" ^ 0.65Mo and A « 0.25, see Wood 1997). 
The amount of dredge-up found in evolutionary calcula¬ 
tions depends sensitively on metallicities, core and enve¬ 
lope masses of the models (Wood 1981). Additionally, as 
pointed out by Frost & Lattanzio (1996), numerical details 
may play an important role. Finally, it is in particular the 
treatment of convection which bears large uncertainties 
concerning the dregde-up efficiency. 

Another problem related to this subject concerns the 
s-process nucleosynthesis. The required neutrons are prob¬ 
ably not released via the ^^Ne(a, n)^®Mg reaction. Instead, 
the s-process is most likely driven by the ^^G(a,n)^®0 
neutron source, raising the question how to mix pro¬ 
tons into carbon-rich layers in order to produce sufficient 
amounts of ^^G . Iben & Renzini (1982) found that the 
protons can diffuse from the bottom of the convective 
envelope into the intershell zone due to semiconvection. 
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However, this scenario is restricted to low metallicities 
(Iben 1983). Often, “standard calculations” artificially in¬ 
gest a given amount of at the onset of the pulse into 
the convectively unstable He shell. However, Straniero et 
al. (1995) have shown that formed during the dredge- 
up phase is burnt under radiative conditions during the 
interpulse phase (for a recent review see Lattanzio 1995). 

These difficulties of standard stellar evolution calcula¬ 
tions to confront the observations often led to the conclu¬ 
sion that mixing outside the formally convective bound¬ 
aries may take place (e.g. Hollowell & Iben 1988, D’Antona 
& Mazzitelli 1996, Wood 1997), or even that only a hy¬ 
drodynamic approach of modelling the H/He interface will 
overcome the drawbacks of the local treatment of convec¬ 
tion (Arlandini et al. 1995). 

In brief the situation can be summarized as follows: 
The mixing length theory (MLT) (Bohm-Vitense, 1958) 
is usually applied to stellar evolution calculations. But 
it does not describe additional (possibly partial) mixing 
beyond the Schwarzschild boundary. To meet the obser¬ 
vations for main sequence stars the instantaneous mixing 
of the convective core is extended beyond the classical 
Schwarzschild boundary by some fraction of the pressure 
scale height. This treatment is commonly referred to as 
overshoot. It has been discussed for example by Along! et 
al. (1993, and ref. therein). The method we present in the 
following has in some aspects similar consequences as for¬ 
mer overshoot treatments but it is based on hydrodynamic 
calculations. 

Freytag et al. (1996) have carried out two-dimensional 
numerical radiation hydrodynamics simulations to study 
the structure and dynamics of a variety of shallow sur¬ 
face convection zones. Their work reveals an independent 
theoretical evidence for the actual existence of extra mix¬ 
ing beyond the boundary of convectively unstable regions. 
The parametrisation of the exponentially declining veloci¬ 
ties of convective elements beyond the classical convective 
border can be applied to one-dimensional stellar evolu¬ 
tion calculations. This treatment leads then to some extra 
partial mixing (Blocker et al. 1997). In the following chap¬ 
ters we confine ourselves to the description of the method 
and the most important consequences. Further implica¬ 
tions and details will be given elsewhere. 


2. The stellar evolution code and the method of 
mixing 

The computations for this study are based on the code de¬ 
scribed by Blocker (1993, 1995) with major modifications. 
Nuclear burning is accounted for via a nucleosynthesis net¬ 
work including 31 isotopes and 74 reactions up to carbon 
burning. We use the most recent opacities (Iglesias et al. 
1992, Alexander & Ferguson 1994). The initial composi¬ 
tion is (Y, Z) = (0.28,0.02), the mixing length parameter 
of the MLT is a = 1.7. For the 3Mq model sequence pre¬ 


sented in this paper a Reimers-type mass loss with 77 = 1 
has been applied. 


2.1. Equation for abundance changes 

The main modification relevant for this study is the in¬ 
troduction of a time-dependent model for mixing in the 
formulation adopted by Langer et al. (1985) for the study 
of semiconvection in massive stars: 


AXi 

fdX,^ 

d 


dt 


^ dMr 

nuc ^ 

[(4,r p) 


with Xi being the mass fraction of the respective isotope. 
The first term on the right-hand side considers the abun¬ 
dance changes due to nuclear burning and the correspond¬ 
ing rates are determined by the nucleosynthesis network. 
The second part describes the mixing of the elements by 
means of a diffusion term given here with respect to the 
mass coordinate M^. The diffusion coefficient D depends 
on the assumed mixing model. For instance, in convective 
unstable regions within the Schwarzschild boundaries D 
is derived from the MLT (see below). 

In our present calculations we make the assumption 
that the burning rates, i.e. {dXi/dt)nuc, do not change 
during one time step. The nuclear network is invoked just 
once for each time step and the resulting burning rates 
enter Eq. 


2.2. The choice of D 

Eq.|| is integrated over the whole star for each isotope. 
For the radiative zones D is obviously Dr = 0. Within the 
convectively unstable regions we follow the prescription of 
Langer et al. (1985): Dc = l/SvJ with I being the mixing 
length and Vc the average velocity of the convective ele¬ 
ments according to the MLT (Bohm-Vitense, 1958). Then 
the diffusion coefficient can be written as 





— g/3(l - /3)Vad(Vrad “ Vad) 

up 



( 2 ) 


where is the pressure scale height, k the opacity, a 
the mixing length parameter, c the velocity of light, p the 
density, g the gravitational acceleration, V^ad and Vad the 
radiative and adiabatic temperature gradient respectively, 
and [3 the gas pressure fraction. 

Finally, for the regions which are immediately adjacent 
to convectively unstable zones we adopt the depth depen¬ 
dent diffusion coefficient derived by Freytag et al. (1996: 
Eq. 9) from their numerical simulations of two-dimensional 
radiation hydrodynamics of time-dependent compressible 
convection : 

Dos = Do exp —Do = vo- Hp, = f ■ Hp, (3) 

tly 


where z denotes the distance from the edge of the con¬ 
vective zone {z = |redge — with r: radius), and is 
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the velocity scale height of the overshooting convective el¬ 
ements at Tedge being proportional to the pressure scale 
height Hp. According to Freytag et al. (1996), Dq has to 
be taken near the edge of the convective zone. Here the ve¬ 
locity field varies only slighty. We take the corresponding 
typical velocity of the convective elements, vg, obtained 
from the MLT in the unstable layers immediately before 
the Schwarzschild border. 





Fig. 1. Diffusion coefficient (Dc, Dos) versus mass coordi¬ 
nate for an TP-AGB model with / = 0.02, 49yr before the 
largest ^He-luminosity (beginning of 12*^ pulse cycle). The 
pulse-driven convective zone (pz) occupies the lower half of 
the intershell zone. The latter extends from 0.587Mq (location 
of the “^He-shell) to 0.607Mq (where H-burning takes place). 
The vertical dotted lines indicate the Schwarzschild boundaries 
(Sb). For the base (Mr = 0.608Mq) of the convective envelope 
(ce) it can not be distinguished from Dos since the H-shell 
is very extended. This can be seen from the dashed-dotted 
line which shows the radius (r). However, this overshoot region 
comprises 5.8 10“®Mq accounted for in 16 mass shells. 

We summarize the choice of D (Fig.^: 

{ Dc instantaneous mixing, ordinary convection 
Dos mixing efficiency declining, overshoot 
region 

Dr no mixing, radiative layers 
2.3. The determination of f 

The parameter / is a measure of the efficiency of the extra 
partial mixing. It defines the velocity scale height of the 
convective elements beyond the boundary of the convec- 
tively unstable zone. As can be seen from Eq. ^the decline 
of Dos is steeper for smaller values of /, thus the larger / 
the further extends the extra partial mixing beyond the 
convective edge. 

Freytag et al. (1996) find from their simulations of 
overadiabatic convective envelopes of A-stars and DA 
white dwarfs / = 0.25 ± 0.05 and 1.0 ± 0.1 respectively. 
These results already reveal that / depends on the stellar 
parameters. For deep envelope (or core) convection / can 
be expected to be considerably smaller since the ratio of 
the Brunt-Vaisala timescales of the stable to the insta¬ 
ble layers decrease with increasing depth, i.e. adiabacity. 


However, hydrodynamical simulations of such deep con¬ 
vection zones are not available yet and we have to rely 
on indirect methods. For example, / ~ 0.02 results in 
the same main-sequence width as calculations by Schaller 
et al. (1992) who applied instantaneous overshoot with a 
parameter of dover/Hp = 0.2. We took this value of / for 
our calculations assuming that it holds also for the deep 
convection zones of AGB stars. Although the width of the 
calculated main sequence depends very sensitively on / 
(A/ = 0.005 corresponds to a displacement of the termi¬ 
nating age main sequence of Alogg ~ 0.08) our experi¬ 
ments with AGB evolution calculations have shown that 
the character of our findings remains qualitatively valid 
even if / is changed by a factor of two. 

3. First results from the calculations 

Applying this treatment of convection to AGB models we 
found the 3’’'^ dredge up for a 3Mq model. The onset was 
already at the 7**^ thermal pulse at Mh = Q.58Mq reach¬ 
ing a dredge-up parameter of A = 0.6 within the next five 
pulses (see Tab. 1). Note, that corresponding calculations 
without extra mixing (standard calculations) did not show 
any dredge up. 

Due to the application of depth-dependent overshoot 
to all convective zones the abundances in the intershell 
zone change in comparison to standard calculations. The 
pulse-driven convective zone extends deeper into the un¬ 
derlying C/O-core already during the first TP and con¬ 
siderable amounts of and are mixed up. Com¬ 
pared to the values given in Tab.|] the abundances of 
the intershell zone in standard calculations is typically 
("'^He/^^C/^®O)=(70/26/l). The comparison of these in¬ 
tershell abundances with observed surface abundances of 
Wolf-Rayet (WR) central stars and PG 1159 stars can be 
used in the future as an additional constraint for /. 


Table 1. Dredge-up data for selected thermal pulses of a 
M = 3Mq sequence: mass of the hydrogen exhausted core 
(Mh), dredge-up parameter (A = AMdup/AMh, AMh: in¬ 
crease of the core mass by nuclear burning during the last pulse 
cycle), total amount of dredged-up material (Mdup), amount 
of dredged up (Mj 3 up(i 2 c))i composition of the intershell 
zone, and content of the ^®C-pocket (Mpsc))- 


TP Nr. 

7 

8 

12 

TJh 

M© 

A 

Mdup 

M© 

*fDUP(12c) 

M© 

(^He/^^C/i^O) 

*f(13c) 

_ 

0.5779 

0.11 

6.3 lO""* 

1.4 >1= lO"'* 
(23/50/25) 
3.9 >1= 10"'^ 

0.5837 

0.12 

7.0 * 10“^ 

1.8 * lO”"* 
(23/49/26) 
3.7* 10~'^ 

0.6073 

0.60 

3.7* 10"® 

1.5* 10"® 
(26/45/27) 
2.0 * 10"'^ 
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Fig. 2. Onset of the 3rd dredge-up during the 12*^ pulse 
cycle. At 0.6072Mq is the bottom of the convective enve¬ 
lope which has already moved into the '*He-dominant layer at 
Mt/Mq ss 0.6071 — 0.6072. The latter is the product of the 
H-shell burning of the previous pulse. This situation is located 
273yr after the beginning of 12**' pulse cycle, i.e. 322yr later 
than shown in Fig. 
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Fig. 3. The convective envelope continues to proceed down¬ 
wards (210yr past Fig.|^ and will shortly hereafter stop. Note 
the different mass range compared to the previous figure. 

This strong increase of in particular causes now a 
corresponding increase of opacity and thereby of the ra¬ 
diative gradient Vrad and leads to semiconvection (Iben 
& Renzini 1982) immediately below the ^He-rich layer. 
In our calculations this shows up as thin convective shells 
which may be connected by their overlaping diffusive tails. 
The effect is that the sharp border between the inter¬ 
shell zone and the '^He-rich layer is wiped out (Fig.||, 
M ss O. 6 O 7 M 0 ). In the following, the envelope convec¬ 
tion continuously progresses downwards, and finally its 


Fig. 4. The downward movement of the convective envelope 
has stopped, the dredge-up is over. A small *®C-pocket has 
formed in the region of previously overlaping diffusion tails of 
th H- and *^C-profile. The snapshot is 17660yr later than the 

Fig. I 

base proceeds through the previously semiconvective zone 
(Fig.|). Typically, 50% of the dredged-up material consist 
of the rest is given by “^He and resp., with 25% 
each (Tab.|^). 

Finally the convective envelope base remains at the 
position of deepest penetration for about fifty years. This 
is the phase where the adjacent H and profiles can 
built up a very thin overlap zone due to the diffusive tail 
of the envelope convection (Uos, Eq.^. At this location in 
the star the *^C-pocket can form (Fig.^. We confirm the 
findings of Straniero et al. (1995) that the built up *^C- 
pocket is burnt under radiative conditions well before the 
next thermal pulse. 

4. Concluding remarks 

We considered overshoot in evolutionary calculations 
based on results of 2-dimensional hydrodynamical simu¬ 
lations of convection. On the AGB, this leads to the for¬ 
mation of a *^C-pocket in the intershell region in contrast 
to calculations based on standard overshoot considera¬ 
tions (i.e. on formal shifts of the Schwarzschild border). 
Additionally, the intershell abundances are considerably 
changed due to extra mixing. As opposed to calculations 
without any overshoot, we find now dregde-up even for 
models with masses as low as SMq. 
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